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INTRODUCTION 

The magnitude of the neural activa- 
tion, and hence the force produced by 
a muscle, depend on the number of 
motor units activated (recruitment) and 
the rates at which motor neurons dis- 
charge action potentials (rate coding). 
Although the recruitment order of motor 
units (size principle) is similar for con- 
tractions during which the force is grad- 
ually increased (ramp contraction) and 
those during which the force is produced 
as fast as possible (see Duchateau and 
Enoka, 2011), rate coding differs between 
the two types of contractions (Desmedt 
and Godaux, 1977a,b; Bawa and Calancie, 
1983). Motor unit discharge rate increases 
progressively during slow ramp contrac- 
tions (Milner-Brown et al., 1973) whereas 
fast contractions involve high instanta- 
neous discharge rate that decreases there- 
after (Desmedt and Godaux, 1977a; Van 
Cutsem et al, 1998). Maximal discharge 
rate during slow isometric ramp contrac- 
tions usually reaches values of 20-50 Hz 
whereas it can attain much higher val- 
ues (>100Hz), albeit briefly, during fast 
contractions (for reviews, see Enoka and 
Fuglevand, 2001; Duchateau and Enoka, 
2011) 

Fast isometric contractions can be per- 
formed in different ways. A first possi- 
bility is to increase force as quickly as 
possible up to a certain level and to 
maintain this force for a few seconds 
(step and hold contraction). An alterna- 
tive way is to produce force as fast as 
possible but to relax the muscle imme- 
diately after the target force is reached. 



Such impulse-like contractions have been 
termed ballistic contractions (Desmedt and 
Godaux, 1977a). Although both contrac- 
tions involved reaching a target force as 
fast as possible, results from our laboratory 
indicate that the maximal rate of torque 
development is ~16% greater for ballistic 
than step and hold contractions (465.2 ± 
17.4 vs. 400.5 ± 20 Nm/s; mean ± SD) 
performed with the ankle dorsiflexor mus- 
cles. Considering the difference in motor 
unit discharge rate between slow and fast 
contractions, these data suggest that ballis- 
tic contractions could be used to assess the 
maximal discharge rate of motor neurons 
in humans. 

MOTOR UNIT DISCHARGE RATE 
DURING BALLISTIC CONTRACTION 

Desmedt and Godaux (1977a) were the 
first to provide a detailed description 
of motor unit discharge in the tibialis 
anterior muscle during ballistic contrac- 
tions. They reported that during bal- 
listic contractions, motor units usually 
began to discharge at high instantaneous 
rates (60-120 Hz) that thereafter declined 
progressively during their successive dis- 
charges, presumably reflecting the ini- 
tial phase of discharge rate adaptation 
observed during repetitive activation of 
motor neurons (Sawczuk et al., 1995; Miles 
et al., 2005). Such discharge pattern has 
been also reported for the first dorsal 
interosseus (Desmedt and Godaux, 1977b) 
and the masseter (Desmedt and Godaux, 
1979), with very brief interspike interval 
(<10ms) mainly observed for the initial 
discharges (Desmedt and Godaux, 1977a; 



Van Cutsem et al, 1998; Van Cutsem and 
Duchateau, 2005). Similar brief interspike 
intervals have also been recorded in the 
flexor carpi radialis during fast and hold 
contractions (Bawa and Calancie, 1983). 
Such high motor unit discharge rates are 
similar to those reported for motor neu- 
rons in animal studies in response to fast 
current injection (Kernell, 1965; Baldissera 
et al, 1987; Sawczuk et al, 1995), and 
should mainly reflect the effect of the 
strong excitatory inputs required to pro- 
duce ballistic contractions. However, these 
very high discharge rates could also be 
influenced by the trajectory of motor 
neuron membrane after de repolariza- 
tion phase (delayed depolarization phase 
and/or after-hyperpolarization period — 
AHP) at the time of the activation (see 
Garland and Griffin, 1999; Kudina and 
Andreeva, 2013). 

TASK-RELATED CHANGES IN 
DISCHARGE RATE 

The discharge characteristics of single 
motor units during ballistic contractions 
can be modulated by the conditions under 
which the action is performed. For exam- 
ple, it has been observed that when a bal- 
listic contraction with the ankle dorsiflex- 
ors was superimposed on a submaximal 
isometric contraction (20-25% of maxi- 
mal force), the average discharge rate for 
the first three interspike intervals was sig- 
nificantly reduced by 22% (89.8 ± 14.6 
vs. 115 ± 20.9 Hz; mean ± SD) com- 
pared with ballistic contractions per- 
formed from a resting state (Van Cutsem 
and Duchateau, 2005). The percentage of 
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motor units that exhibited discharges rate 
above 200 Hz at the onset of the activa- 
tion was also diminished (6.2 vs. 15.5%). 
Interestingly, the instantaneous discharge 
for the first interspike interval was much 
reduced (-37%) than the second (-18%) 
and third (—8%) intervals. This lower 
motor unit discharge rate during super- 
imposed ballistic contractions was accom- 
panied by a decrease in the maximal 
rate of force development (~16%). The 
slower rate of force development and 
reduced motor unit discharge rate dur- 
ing the superimposed ballistic contrac- 
tions are, however, abolished when a brief 
silent period (usually called "premotor 
silent period") was observed at the transi- 
tion between the pre-activation (sustained 
contraction) and ballistic actions (Van 
Cutsem and Duchateau, 2005). A simi- 
lar observation has been reported when a 
brief voluntary agonist relaxation (deac- 
tivation) was inserted between the sus- 
tained and the ballistic action (Duchateau 
and Baudry, 2012). These silent periods 
(unintentional and voluntary) are thought 
to enable motor neurons to achieve a 
non-refractory state leading to a more 
synchronous recruitment and a greater 
discharge rate of motor units during 
the subsequent ballistic action (Tsukahara 
et al, 1995; Van Cutsem and Duchateau, 
2005). The changes in maximal discharge 
rate achieved during ballistic contractions 
with initial conditions likely reflect the 
history-dependent changes of motor neu- 
ron excitability (Heckman and Enoka, 
2012), and on a functional point of view 
supports the association between the max- 
imal motor unit discharge rate and the rate 
of force development. 

LONG-TERM CHANGES IN DISCHARGE 
RATE 

A way to further investigate this asso- 
ciation consists of studying long-term 
changes in the maximal discharge rate of 
human motor units, such as those occur- 
ring in response to training and ageing. 
For example, Van Cutsem et al. (1998) 
reported that 3 months of ballistic con- 
tractions of the ankle dorsiflexor muscles 
against a moderate load (30-40% MVC) 
enhanced the maximal rate of force devel- 
opment by 82% during ballistic contrac- 
tions. Although no change was observed 
in the recruitment order of motor units, 



the average discharge rate of the first 
four action potentials increased by 38% 
after training (96.3 ± 39.5 vs. 69.9 ± 
30.8 Hz; mean ± SD). The increase in 
discharge rate was significantly less for 
the first (+86%) and second (+70%) 
than the third (+124%) interspike inter- 
vals. In addition, training increased the 
number of motor units (5-33%) exhibit- 
ing discharges above 200 Hz at the onset 
of activation. Because the average time 
to peak force of motor unit mechanical 
responses was not significantly modified, 
the increase in the rate of force devel- 
opment during the ballistic contractions 
was mainly due to adaptation in motor 
unit discharge rate. Potential mechanisms 
that may explain the changes in motor 
unit discharge rate should involve dif- 
ferent loci along the corticospinal path- 
way. Although some of these changes can 
occur at supraspinal level (Schubert et al., 
2008), part of the adaptations presumably 
involve changes in the intrinsic proper- 
ties of motor neurons, as observed after 
endurance training in rats (Gardiner et al., 
2006). 

In contrast to training, the ageing pro- 
cess induces a decline in the speed-related 
capacity of individuals. For example, the 
maximal rate of force development during 
ballistic contractions performed with the 



ankle dorsiflexor muscles was significantly 
lower by 48% in elderly (71-84 year) 
than in young adults (~20 year) (Klass 
et al., 2008). This age-related change was 
accompanied by a clear decline in the 
average motor unit discharge rate. As the 
decrease was less pronounced for the first 
(-19%) than for the second (-28%) 
and third (—34%) interspike intervals, 
this means that the aged motor units 
cannot sustain a high discharge rate dur- 
ing successive discharges. In addition, the 
percentage of motor units that exhib- 
ited initial discharges above 200 Hz was 
reduced (—45%) in elderly compared with 
young adults. As the rate of force devel- 
opment during electrically evoked con- 
tractions, that by-pass motor neurons 
activation, is less reduced than those 
during ballistic voluntary contractions, the 
decline in maximal motor unit discharge 
rate should significantly contribute to limit 
the performance of fast voluntary contrac- 
tions with ageing. The age-related prolon- 
gation in the duration of motor neuron 
after hyperpolarization, as observed in 
the human biceps brachii by Piotrkiewicz 
et al. (2007), could be a relevant candi- 
date to explain, at least in part, the reduced 
maximal rate of motor unit discharge 
during ballistic contractions in elderly 
adults. 
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FIGURE 1 | Simulation of the relation between motor unit discharge rate and maximal rate of 


force development for the 1st, 100th, and 200th motor unit (MU) of a pool of 200 units in the 


tibialis anterior muscle. The simulation was based on a model developed by Fuglevand et al. 


(1993) with the inclusion of the spike-triggered average forces for motor units published by Van 


Cutsem et al. (1998). The force 


generated by each motor unit was simulated for 4 successive 


discharges generated at constant frequencies ranging from 10 to 500 Hz before the first derivative 


was computed to obtain the maximal rate of force development. 
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MODELING THE RELATION BETWEEN 
MOTOR UNIT DISCHARGE RATE AND 
RATE OF FORCE DEVELOPMENT 

To further analyse the effect of a change 
in discharge rate on the maximal rate of 
force development, isometric force pro- 
duced by single motor units was simu- 
lated from a model that contains a pool 
of 200 units (Fuglevand et al., 1993; 
Duchateau and Enoka, 2002). To that pur- 
pose, mechanical properties (peak force 
and time to peak force) of motor units 
obtained from the spike-triggered aver- 
aging method in the tibialis anterior 
(Van Cutsem et al., 1998) were inserted 
into the model. Data indicated that an 
increase in discharge rate up to 100- 
200 Hz augmented substantially the rate 
of force development for all units of 
the pool (Figure 1). Nonetheless, further 
increase in discharge rate has less influ- 
ence excepted for the faster units (MU 
100 and MU 200) of the pool, reflect- 
ing difference in speed-related properties 
between low- and high threshold motor 
units. These simulated data underscore 
the critical role of maximal motor unit 
discharge rate on the ability to rapidly 
develop force. 

CONCLUDING REMARKS 

Together, experimental and simulated data 
indicate that a high initial motor unit dis- 
charge rate at the onset of a fast contrac- 
tion plays a critical role to reach a high 
rate of force development. Furthermore, 
and because the instantaneous discharge 
rates of motor units at the onset of ballistic 
contractions are much greater than those 
recorded during slow contractions and 
not yet influenced by history-dependent 
effects, ballistic contractions from a rest- 
ing state can be used to assess the maximal 
motor neuron discharge rate in human. 
Nonetheless, as the acquisition of a simple 
motor task such as index finger abduc- 
tion requires up to ~300 repetitions to 
reach maximal acceleration capability (Lee 
et al., 2010), subjects must be familiarized 
beforehand with ballistic contractions of 
the muscle under study. 
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